Maria Figures Final
Results from Marine Heatwave Simulations 2023
Loading the Data
Data has been freshly updated in Aufust 2023 to be consistent with an error/bug found in the SST data.
Locations to Files: Here is the link to the resources used to create these figures:
Species Data Key
Dominant Species - Data Wrangle:
Couple data wrangling steps here: - make row and column identifiers for what the season was, and what the year was that matches the heatwave trial - add one year of control as a lead-in for plots - join the control biomass as a column so we can caluclate difference and percent difference from what it would be without a heatwave disturbance - do some benchmarking against the taxa specific intra-annual fluctuation in biomass
Prepare the absolute biomass data for the dominant species:
Dominant Species - Seasonal Cycle Control
NOTES: Do a similar cleanup to the heatwave results, with some minor tweaks want all the columns to match so they can append. year = 0, hw_season is either none, or we can repeat it 4 times to filter with the data or we can just plop it on top whenever we want it later.
Dominant Species Data Wrangling:
1. Isolate the HW Season Responses:
2. Lead with 1-Year of Control:
3. Calculate Difference from Control:
Fig. 1: Abiotic Conditions
The following figures show the annual cycles of abiotoc conditions as their own plots. This is the context for how each seasonal cycle progresses within the year, and what the plankton community are responding to.
Fig. 2 Control Conditions
Make the control cycles match the design of Marias python figures X axis: time (daily for a, seasonal for b and c). Y axis: absolute biomass. In subplots (a): dashed line shows the end of each season. In subplots (c): Bars: temperature optima. Stacked bars: absolute biomass concentrations of the different size groups (color bar for the size range of groups from small to big. You can use another colorbar 😊).
Subtitle for each season on x axis Top to bottom: active, passive, protists
Fig. 2 Complete
Fig 3. Size Bin Response
Absolute Biomass concentration anomalies (HW/afterHW – preHW) for the biomass size bins of (a) protists, (b) passive and (c) active copepod feeders for the winter, autumn, spring and summer HW scenarios.
Fig 3 Option 2.
These focus in on one specific seasonal heatwave. This makes it easier to see progression within one specific season at the expense of being able to see the annual cycles
Fig 4. Dominant Species Biomass
Absolute seasonal biomass concentration anomalies of the dominant groups of (a) protists, (b) passive and (a) active copepod feeders for the autumn HW scenario during the heatwave and 9 years after the heatwave.
Fig 5. Shannon Diversity
Annual Shannon Diversity Index for (a) protists, (b) passive and (c) active copepod feeders for the winter, spring, autumn and winter HW scenarios. y axis: years. Year 0: pre-heatwave, Year 1: heatwave (in red) Years 2-9: after-heatwave.
Note for Adam: Here we will have the figure you have under the Annual Change within a Season section. Same ranges for the Shannon index for all groups (2.4- 6.5). X axis: can you please change it to 0 1 2 3 4 5 6 7 8 9 than 0 2 4 6 8 years?
Fig 6. Biomass Anomalies for Differential Exposure
Comment:
Take just the values starting with: during the heatwave, and then 7 seasons after (2 years total)
Biomass anomalies (HW protists or copepods only - HW on both protists and copepods) for all seasons for the seasonal HWs and the 7 seasons afterwards. X axis: size bins, y axis: absolute biomass anomalies. The median value is presented with a square and the Standard deviation with the bar line.
Comments
Figure 4 needs to show relative dominance. Not the change in biomass from the control. Needs to highlight how the top x species within the different groups are competing/jockeying Actives: top 5 Passives: top 5 Protists: top 10 or 5
Figure 6, filter out the season before the heatwave, and then take the heatwave season and the following 7 seasons after that to show the data.